Safe Cooperative Robot Dynamics on Graphs* 



Robert W. Ghrist 1 " Daniel E. Koditschek* 

February 1, 2008 



Abstract 

This paper initiates the use of vector fields to design, optimize, and implement reactive schedules 
for safe cooperative robot patterns on planar graphs. We consider Automated Guided Vehicles (AGV's) 
operating upon a predefined network of pathways. In contrast to the case of locally Euclidean configu- 
ration spaces, regularization of collisions is no longer a local procedure, and issues concerning the global 
topology of configuration spaces must be addressed. The focus of the present inquiry is the achievement 
of safe, efficient, cooperative patterns in the simplest nontrivial example (a pair of robots on a Y-network) 
by means of a state-event heirarchical controller. 



1 Introduction 

Recent literature suggests the growing awareness of a need for "reactive" scheduling wherein one desires not 
merely a single deployment of resources but a plan for successive re-deployments against a changing environ- 
ment . But scheduling problems have been traditionally solved by appeal to a discrete representation of 
the domain at hand. Thus the need for "tracking" changing goals introduces a conceptual dilemma: there is 
no obvious topology by which proximity to the target of a given deployment can be measured. In contrast 
to problems entailing the management of information alone, problems in many robotics and automation set- 
tings involve the management of work — the exchange of energy in the presence of geometric constraints. 
In these settings, it may be desirable to postpone the imposition of a discrete representation long enough to 
gain the benefit of the natural topology that accompanies the original domain. 

This paper explores the use of vector fields for reactive scheduling of safe cooperative robot patterns on 
graphs. The word "safe" means that obstacles — designated illegal portions of the configuration space — are 
avoided. The word "cooperative" connotes situations wherein physically distributed agents are collectively 
responsible for executing the schedule. The word "pattern" refers to tasks that cannot be encoded simply in 
terms of a point goal in the configuration space. The word "reactive" will be interpreted as requiring that 
the desired pattern reject perturbations: conditions close but slightly removed from those desired remain 
close and, indeed, converge toward the exactly desired pattern. 



1.1 Setting: AGV's on a Guidepath Network of Wires 

An automated guided vehicle (AGV) is an unmanned powered cart "capable of following an external guidance 
signal to deliver a unit load from destination to destination" where, in most common applications, the 
guidepath signal is buried in the floor ||. Thus, the AGV's workspace is a network of wires — a graph. The 
motivation to choose AGV based materials handling systems over more conventional fixed conveyors rests not 
simply in their ease of reconfigurability but in the potential they offer for graceful response to perturbations in 
normal plant operation. In real production facilities, the flow of work in process fluctuates constantly in the 
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face of unanticipated workstation downtime, variations in process rate, and, indeed, variations in materials 
transport and delivery rates |7]]. Of course, realizing their potential robustness against these fluctuations in 
work flow remains an only partially fulfilled goal of contemporary AGV systems. 

Choreographing the interacting routes of multiple AGVs in a non-conflicting manner presents a novel, 
complicated, and necessarily on-line planning problem. Nominal routes might be designed offline but they 
can never truly be traversed with the nominal timing, for all the reasons described above. Even under normal 
operating conditions, no single nominal schedule can suffice to coordinate the workflow as the production 
volume or product mix changes over time: new vehicles need to be added or deleted and the routing scheme 
adapted. In any case, abnormal conditions — unscheduled process down times; blocked work stations; failed 
vehicles — continually arise, demanding altered routes. 

The traffic control schemes deployed in contemporary AGV systems are designed to simplify the real-time 
route planning and adaptation process by "blocking zone control" strategies. The workspace is partitioned 
into a small number of cells and, regardless of the details of their source and destination tasks, no two AGVs 
are ever allowed into the same cell at the same time [|| . Clearly, this simplification results in significant loss 
of a network's traffic capacity. 

In this paper, we will consider a centralized approach that employs dynamical systems theory to focus 
on real-time responsiveness and efficiency as opposed to computational complexity or average throughput. 
No doubt, beyond a certain maximum number of vehicles, the necessity to compute in the high dimensional 
configuration space will limit the applicability of any algorithms that arise. However, this point of view 
seems not to have been carefully explored in the literature. Indeed, we will sketch some ideas about how 
an approach that starts from the coupled version of the problem may lend sufficient insight to move back 
and forth between the individuals' and the group's configuration spaces even in real time. For the sake 
of concreteness we will work in the so-called "pickup and delivery" (as opposed to the "stop and go" 
paradigm of assembly or fabrication, and we will not be concerned with warehousing style AGV applications. 

1.2 Contributions of the Paper 

The paper is organized as follows. In we review fundamental facts about the topology of graphs, with 
which we define the class of edge point fields — locally defined dynamics that realize single letter patterns. 
These act collectively as a toolbox from which to build a hybrid controller for achieving arbitrary patterns 
with a single AGV. This represents a slight generalization of the scheme the second author and colleagues 
have proposed in 

The problem of dynamics and control on non-trivial graphs is then considered in beginning with a 
detailed discussion of a natural intrinsic coordinate system in which to frame the configuration space. We 
present a topological analysis of the configuration space for a pair of AGVs on a Y-shaped graph — the 
simplest nontrivial situation. Here, a clarification of the configuration space presentation leads easily to a 
vector field construction that brings all initial conditions of two robots on the graph to any desired pair of 
goal points while guaranteeing safety (i.e., no collisions along the way). The desire for a more decoupled 
controller — the hope of an "interleaving" of otherwise independent individual patterns — impels a revised 
approach to safe navigation leading to the construction of a vector field that enables the AGVs to "dance" 
about one other at a vertex. 

The dynamical features of this circulating field are suggestive of a hybrid construction that would allow 
multiple independent patterns to be safely interleaved. We proceed in §^ to construct a I2-symbol grammar 
of so-called "monotone" cycles: those patterns which exclude multiple robots on a single edge of the graph. 
The goal of this excursion is to tune limit cycles visiting various docking stations in such a way as to be 
optimal with respect to certain notions of distance or performance. 

We complete our treatment of this fundamental example by synthesizing the grammar results into a 
state-event controlled hybrid system for achieving cooperative patterns. Appendix A is included to place on 
a rigorous foundation the use of vector fields on graphs and configuration spaces thereof. 



2 Notation and Background 
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2.1 Graph Topology 

A graph, T, consists of a finite collection of O-dimensional vertices V:—{vi}^, and 1-dimensional edges 
£:={ej}f I assembled as follows. Each edge is homeomorphic to the closed interval [0, 1] attached to V along its 
boundary points {0} and We place upon V the quotient topology given by the endpoint identifications: 

Neighborhoods of a point in the interior of are homeomorphic images of interval neighborhoods of the 
corresponding point in [0, 1], and neighborhoods of a vertex Vi consist of the union of homeomorphic images 
of half-open neighborhoods of the endpoints for all incident edges. 

The configuration spaces we consider in §|| and following are self- products of graphs. The topology of 
r x T is easily understood in terms of the topology of T as follows [0. Let (x, y) G T x T denote an ordered 
pair in the product. Then any small neighborhood of (x, y) within T x T is the union of neighborhoods 
of the form M(u) x J\f(v), where Af(-) denotes neighborhood within T. In other words, the products of 
neighborhoods form a basis of neighborhoods in the product space. 

Given a graph, T, outfitted with a finite number N of non-colliding AGV's constrained to move on T, 
the (labeled) configuration space of safe motions is defined as 

C:=(Tx ... xT)-Af(A), (1) 

where A:={(xi) G V for some j ^ k} denotes the pairwise diagonal and Af(-) denotes 

(small) neighborhood. 

For general graphs, the topological features of C can be extremely complicated. We do not treat the 
general aspects of this problem comprehensively in this paper; rather, we restrict attention to the simplest 
nontrivial example which illustrates nicely the relevant features present in the more general situation. The 
topological characteristics of general configuration spaces on graphs is treated in ||, [| . Mathematically, it is 
usually most interesting to pass to the quotient of C by the action of the permutation group on N elements, 
thus forgetting the identities of the AGV elements; however, as such spaces are almost completely divorced 
from any applications involving coordinated transport, we work on the "full" configuration space C. 

In order to proceed, it is necessary to clarify what we mean by a vector field on a simplicial complex that 
fails to be a manifold. This is a nontrivial issue: for example, in the case of a graph, the tangent space to a 
vertex with incidence number greater than two is not well-defined. We defer a more detailed discussion of 
these statements to Appendix A. The essential difference is that we construct semiflows: flows which possess 
unique forward orbits. 

2.2 Edge Point Fields for Single AGV Control 

In the context of describing and executing patterns or periodic motions on a graph, one desires a set of 
building blocks for moving from one goal to the next. We introduce the terminology and philosophy for 
constructing patterns by way of the simplest possible examples: a single AGV on a graph. This avoids the 
additional topological complications present in the context of cooperative motion. 

We thus introduce the class of edge point fields as a dynamical toolbox for a hybrid controller. Given a 
specified goal point g G ej within an edge of T, an edge point field is a locally defined vector field X g on T 
with the following properties: 

Locally Defined: X g is defined on a neighborhood N(ej) of the goal-edge ej within the graph topology. 
Furthermore, forward orbits under X g are uniquely defined. 

Point Attractor: every forward orbit of X g asymptotically approaches the unique fixed point g G ej^ 

Navigation-Like: X g admits a C° Lyapunov function, <E> 9 : T — > R. 

The following existence lemma (whose trivial proof we omit) holds. 

1 We will assume away in the sequel the possibility of "homoclinic" edges whose boundary points are attached to the same 
vertex. 

2 When it is not clear from the context, we shall denote the goal point achieved by an edge point flow as g(X g ) = {g}. 
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Lemma 1 Given any edge ej C T which is contractible within T, there exists an edge point field X g for any 
desired goal g £ ej. 

The only occasion for which an edge ej is not contractible in T is in the "homoclinic case" when both 
endpoints of ej are attached to the same vertex, forming a loop. In such instances, one may avoid the 
problem by subdividing the edge to include more vertices, which is very natural in the setting of this paper, 
since vertices correspond to workstations along a path. 



2.3 Discrete Regulation of Patterns 



We adopt the standard framework of symbolic dynamics 1 13 . By an excursion on a graph is meant a (possibly 
infinite) sequence of edges from the graph, E = a x . . . ei N ... 6 £ z , having the property that each pair of 
contiguous edges, and ei j+1 share a vertex in common. The set of excursions forms a language, C: the 
so-called subshift on the alphabet defined by the named edges (we assume each name is unique) p3| . The 
shift operator, a, defines a discrete dynamical system on the set of excursions, mapping the set of infinite 
sequences into itself by decrementing the time index. An M -block extension of the original language arises 
in the obvious way from grouping together each successive block of M contiguous letters from an original 
sequence, and it is clear how a induces a shift operator, a M on this derived set of sequences. 

Given a legal block, B = . . . ei M € £, we will say that an excursion realizes that pattern if its M-block 
extension eventually reaches the "goal" BBBBB . . . under the iterates of a M . In other words, after some 
finite number of applications of er, the excursion consists of repetitions of the block B (terminating possibly 
with the empty edge). 

In a previous paper Q , the second author and colleagues introduced a very simple but effective discrete 
event controller for regulating patterns on graphs from all reachable initial edges by pruning the graph back 
to a tree (imposing an ordering). Of course, this simple idea has a much longer history. In robotics it was 
introduced in JT4[ as "pre-image backchaining;" pursued in [|l5| as a method for building verifiable hardened 
automation via the metaphor of a funneling; and in Q as a means of prescribing sensor specifications from 
goals and action sets. In the discrete event systems literature an optimal version of this procedure has been 
introduced in and a generalization recently has been proposed in | fL8| . 

Let £°:=B c £ denote the edges of T that appear in the block of letters specifying the desired pattern. 
Denote by 

£ n+1 C £ - [j£ k 

k<n 

those edges that share a vertex with an edge in £ n but are not in any of the previously defined subsets. This 
yields a finite partition of £ into "levels," {£ p }p = o, such that for each edge, ef S £ p , there can be found 
a legal successor edge, e v ~ x e £ p ~ 1 , such that e\e p ~ x £ £ is a legal block in the language. Note that we 
have implicitly assumed £° is reachable from the entire graph — otherwise, there will be some "leftover" 
component of £ forming the last cell in the partition starting within which it is not possible to achieve the 
pattern. Note as well that we impose some ordering of each cell £ p = {ef}*^: the edges of £ = B are 
ordered by their appearance in the block; the ordering of edges in higher level cells is arbitrary. 

We may now define a "graph controller" law, G:£ — > £ as follows. From the nature of the partition {£ p } 
above, it is clear that the least legal successor function, 

J i + 1 mod M : p = 



min{j < M p : e\e p r l G £} :p>0 
is well-defined. From this, we construct the graph controller: 

G(eD:=4, b V (3) 

It follows almost directly from the definition of this function that its successive application to any edge leads 
eventually to a repetition of the desired pattern: 



Proposition 2 The iterates of G on £ achieve the pattern B. 
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2.4 Hybrid Edge Point Fields 

A semiflow, (X) 1 , on the graph induces excursions in C parametrized by an initial condition as follows. The 
first letter corresponds to the edge in which the initial condition is located (initial conditions at vertices are 
assigned to the incident edge along which the semiflow points). The next letter is added to the sequence by 
motion through a vertex from one edge to the next. 

We will say of two edge point fields, X\, X2 on a graph, T, that X\ prepares X2, denoted X\ y X2, if the 
goal of the first is in the domain of attraction of the second, g(Xi) C A/^X). Given any finite collection of 
edge point fields on T, we will choose some < a < 1 and assume that their associated Lyapunov functions 
have been scaled in such a fashion that X\ >~ X2, implies ($1) 1 [0, a] C A/^X). In other words, an a 
crossing of the trajectory $1 o (Xif signals arrival in A^X). 

Suppose now that for every edge in some pattern block, G £°, there has been designated a goal point, 
<7°, along with an edge point field X® taking that goal, g{X®) = g®. Assume as well that the edge point field 
associated with each previous edge in the pattern prepares the flow associated with the next edge, in other 
words, using the successor function (|^) we have, 

g(X°)cAf(xl 



r 4i) 

Now construct edge point fields on all the edges of T such that the tree representation of their >- relations 
is exactly the tree pruned from the original graph above — namely we have 



g(Aj)cA^) 



We are finally in a position to construct a hybrid semi-flow on T. This feedback controller will run the 
piece-wise smooth vector field, x — X, as follows 

fX? :x G and > a 

It is clear from the construction that progress from edge to edge of the state of this flow echoes the graph 
transition rule G, constructed above. 

Proposition 3 The edge transitions induced by the hybrid controller (Q) are precisely the iterates of the 
graph map, G, ^) in the language, C. 



3 The Y-Graph 

We now turn our attention to the safe control of multiple AGV's on a graph work-space via vector fields. 
For the remainder of this work, we consider the simplest example of a non-trivial configuration space: that 
associated to the Y-graph, T, having four vertices {vi}l and three edges {e^ }f , as illustrated in Figure [I]. 
Although this is a simple scenario compared to what one finds in a typical setting, there are several reasons 
why this example is in many respects canonical. 

1. Simplicity: Any graph may be constructed by gluing X-prong graphs together for various K. The 
K = 3 model we consider is the simplest nontrivial case and is instructive for understanding the 
richness and challenges of local cooperative dynamics on graphs. 

2. Genericity with respect to graphs: Graphs which consist of copies of T glued together, the trivalent 
graphs, are generic in the sense that any nontrivial graph may be perturbed in a neighborhood of the 
vertex set so as to be trivalent. For example, the 4-valent graph resembling the letter 'X' may be 
perturbed slightly to resemble the letter 'FT — a trivalent graph. An induction argument shows that 
this is true for all graphs. Hence, the dynamics on an arbitrary graph are approximated by patching 
together dynamics on copies of T. 
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3. Genericity with respect to local dynamics: Finally, pairwise local AGV interactions on an ar- 
bitrary graph restrict precisely to the dynamics of two agents on T as follows. Given a vertex v of a 
graph r, assume that two AGV's, x and y, are on different edges e\ and ei incident to v and moving 
towards v with the goal of switching positions. A collision is imminent unless one AGV "moves out of 
the way" onto some other edge e% incident to v. The local interactions thus restrict to dynamics of a 
pair of AGV's on the subgraph defined by {u; ei, e%, e^}. Hence, the case we treat in this paper is the 
generic scenario for the local resolution of collision singularities in cooperative dynamics on graphs. 

3.1 Intrinsic Coordinates 

The configuration space of N points on T is a subset of the iV-fold cartesian product T x T x • • • x T. 
Since each graph which is physically relevant to the setting of this paper is embedded in a factory floor 
or ceiling and thus planar, the configuration space C N (T) embeds naturally in R 2N . We wish to modify 
this embedding to facilitate both analysis on and visualization of the configuration space. We will present 
alternate embeddings in both higher and lower dimensional Euclidean spaces for these purposes. 

We begin with representing the configuration space within a higher-dimensional Euclidean space via a 
coordinate system which is intrinsic: it is independent of how the graph is embedded in space. We illustrate 
this coordinate system with the Y-graph T, noting that a few simple modifications yields a coordinate scheme 
for more general graphs. 

Let {ei]\ denote the three edges in T, parametrized so that is identified with [0, 1] with each {0} at 
the center vq of T. Any point in T is thus given by a vector x in the {e^} basis whose magnitude |x| G [0, 1] 
determines the position of the point in the direction. Denote by i(x) the value of i so that x — \x\e L ^ x y 
This parametrization embeds T as the positive unit axis frame in R 3 . Likewise, a point in C is given as a pair 
of distinct vectors (x, y), i.e., as the positive unit axis frame in R 3 cross itself sitting inside of R 3 x R 3 = R 6 . 
We have thus embedded the configuration space of two distinct points on T in the positive orthant of R e . It 
is clear that one can embed the more general configuration space of N points on T in R 3N in this manner. 

This coordinate system is particularly well-suited to describing vector fields on C and in implementing 
numerical simulations of dynamics, as the coordinates explicitly keep track of the physical position of each 
point on the graph. 

3.2 A Topological Analysis 

More useful for visualization purposes, however, is the following construction which embeds C within R 3 . 




Figure 1: [left] The Y-graph T; [right] the configuration space C embedded in R 3 . 



Theorem 1 The configuration space C associated to a pair of AGV's restricted to the Y-graph T is home- 
omorphic to a punctured disc with six 2-simplices attached as in Figure [3. 
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Proof: Recall that C consists of pairs of distinct vectors (x, y) in intrinsic coordinates. Restrict attention 
to the subspace DcC defined by 

V:={(x,y)eC:i(x)^i(y)}, (5) 

where an undefined index is considered to be different than one which is defined. Thus, T> consists of those 
configurations for which both AGV's do not occupy the same edge interior. 

The set V has a natural cellular decomposition as follows. There are 2 AGV's and 3 edges in T; hence, 
there are 3-2 = 6 cells T>ij C T> where i := i(x) ^ i(y) =: j. Since (the closure of) each edge in T is 
homeomorphic to [0, f] (determined by | • |), the cell is homeomorphic to ([0, f] x [0, f]) — {(0, 0)}, where, 
of course, the origin (0, 0) is removed as it belongs to the diagonal A. A path in T> can move from cell to 
cell only along the subsets where the index of one AGV changes: e.g., \x\ = or \y\ = 0. Thus, the edges 
{0} x (0, 1] and (0, f] x {0} of the punctured square T>ij are attached respectively to T>kj and X>i,fc, where 
k is the unique index not equal to i or j. 

Furthermore, each 2-cell T>ij has a product structure as follows: decompose T>ij along lines of constant 

9 := tan -1 fjfjj- It i s dear that 9 is the angle in the unit first quadrant in which T>i j sits. Hence, each 

T>i t j is decomposed into a product of a closed interval Sij := 9 6 [0, n/2] (an 'angular' coordinate) with the 
half-open interval (0, 1] (a 'radial' coordinate). As this product decomposition is respected along the gluing 
edges, we have a decomposition of all of V into the product of (0, 1] x S, where 5* is a cellular complex 
given by gluing the six segments $y end-to-end cyclically along their endpoints. The set S is a f -manifold 
without boundary since each S% t j is a closed interval, each of whose endpoints is glued to precisely one other 
Si.j. Hence, by the classification of f-manifolds, S is homeomorphic to a circle. We have thus decomposed 
T> as the cross product of a circle with (0, 1] — a punctured unit disc. 

The complement of V in C consists of those regions where l(x) — i(y). For each i — L.3, the subset 
of C where l(x) — i(y) — i is homeomorphic to ((0, f] x (0, 1]) — {\x\ = \y\}: this consists of two disjoint 
triangular "fins." A total of six such fins are thus attached to T> along the six edges where |a;| or \y\ = 0. In 
the coordinates of the product decomposition for T>, these fins emanate along the radial lines where 9 equals 
zero or 7r/2, yielding the topological space illustrated in Figure [|. o 

The precise drawing of Figure [l] represents this punctured disc T> as a hexagon-shaped complex with a 
punctured center: this follows from the cellular structure of T> as being built from six squares sewn together. 

Corollary 4 Given any pair of goals g'=(gi, 32) where g\ and g2 live on different branches of T, there 
exists a navigation function (of class piecewise real-analytic) generating a semifiow which sends all but a 
measure-zero set of initial conditions to g under the gradient semifiow. 

Proof: The subset V C C is homeomorphic to a punctured disc S x (0, 1] , and may be easily compactified 
to an annulus 5* x [0, 1] by removing an open neighborhood of the diagonal. Then, the conditions for the 
theorems of Koditschek and Rimon |fl2f are met, since an annulus is a sphereworld. Hence, not only does a 
navigation function $ on this subspace exist, but an explicit procedure for determining $ is given [pL2|| . One 
may then extend $ to the remainder of C as follows: choose a point (cc, y) on the fin and define 

$(x,y):={ , (6) 

so that $ increases sharply along the fins. This directs the gradient flow to monotonically "descend" away 
from the diagonal and onto T>. Note that T> is forward-invariant under the dynamics, and, that upon 
prescribing the vector field on the fins to point transversally into T>, we have defined a semifiow, and hence 
a well-defined navigational procedure. o 

This result is very satisfying in the sense that it guarantees a navigation function by applying existing 
theory to a situation which, from the definition alone, would not appear to be remotely related to a sphere- 
world. However, a deeper analysis of configuration spaces of graphs Q reveals that for more than two AGV's, 
the configuration space is never a sphereworld.^ Hence, we consider an alternate solution to the problem of 

3 The configuration space of a graph turns out to be aspherical: there are no essential closed spheres of dimension larger 
than one. 
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realizing compatible goals by means of a vector field on the configuration space. This method is adaptable 
to more complicated settings. 

3.3 Example: a Circulating Flow 

Before proceeding with a general scheme for controlling two agents on the Y-graph, we present a simple 
example of a vector field on the configuration space which can be used to regularize collisions about a 
generic trivalent vertex. Theorem ^ and Figure [l] suggest a natural circulating flow on the configuration 
space C which has the effect of inducing a "dance" between the pair of AGV's which cycles through all 
combinations of distinct point goals. 

Theorem 2 There exists a piecewise- smooth vector field X on C which has the following properties: 

1. X defines a nonsingular semiflow on C; 

2. The diagonal A is repelling with respect to X ; and 

3. Every orbit of X approaches a unique attracting limit cycle on C which cycles through all possible 
ordered pairs of distinct edge-states. 

Proof: Denote by V that portion of the configuration space which corresponds to the AGV's being 
on distinct edges of the graph: as proven earlier, T> is homeomorphic to a punctured disc. The intrinsic 
coordinates on the configuration space C is illustrated in Figure [^, where only T> is shown for illustration 
purposes. The reader should think of this as a collection of six square coordinate planes, attached together 
pairwise along axes with the origin removed (this is actually an isometry for the natural product metric). 
The six triangular fins are then attached as per Figure El. 




Figure 2: [left] The coordinate system on the unfinned region V of C; [right] The circulating flow with a 
typical orbit. 

Recal that any point in the graph is represented as a vector x = \x\ei for some i. Denote by the unit 
tangent vector in each tangent space T x ei pointing in the positive (outward) direction towards the endpoint 
Vi. The vector field we propose is the following: given (x, y) E C, 



1. If l(x) = i(y) then 

J ± = ~W^) \ < |x| < |«| 



(7) 

x = \x\(l — \x\)e,( x ) 1 _ _ | . _ . | 

• _ V r <\v\ < M 



V = - \x\e L ( y) 
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2. If l{x) = i{y) + 1 or |x| = then 



X = Z\ l{V) t^ 0<\x\<\y\ 

V = 1 2/1 ( 1 - \y\)e L (y) f 

x = 1x1(1 — b|)e,(' x i | _ | . . , | 

V ? < \y \ < \x\ 

y = -fK(i/) J 



(8) 



3. If t(y) = i(x) + 1 or \y\ = then 



x = |x|(l - laiDe,,^) 
2) = Me( t (x)+i) 



0<|x[ < |y| 

(9) 

0< \y\ < \x\ 



Note that all addition operations on i{x) and i{y) are performed mod three. 

The vector field is nonsingular as follows: if |x||y| ^ 0, then the vector field is by inspection nonsingular. 
If |x| = 0, then \y\ > = \x\ since the points are distinct. It then follows from Equation (||) that the vector 
field on this region has d\x\/dt = \y\ ^ 0. A similar argument holds for the case where \y\ = 0. 

The vector field defines a semiflow as follows: on those regions where ^ \x\ ^ \y\ ^ 0, the vector field 
is smooth and hence defines a true flow. Along the lines where |x| = \y\, the vector field is only C°, but 
nevertheless is constructed so as to define unique solution curves; hence the region V, where t(x) ^ i(y), 
is invariant under the flow. Finally along the branch line curves where |x| = or \y\ = 0, the vector field 
points into the the branch lines from the fins, implying that the dynamics is a semiflow (see the remarks in 
Appendix A). 

This vector field admits a C° Lyapunov function $ : C — > [0, 1) of the form 

\ 1 - max{|x|, |y|} : t(x) ^ i{y) K ' 

From Equation (0), one computes that on the fins (l(x) = o(y)), 

<0, (11) 



~dl 



d\x\ _ d\y\_ 
[ dt dt 1 



since here \x\ ^ \y\. Furthermore, on the disc V (i(x) ^ t(y)), changes as 

§ = *(*-!). (12) 

Hence, $ strictly decreases off of the boundary of the disc 

dV:= {(x, y) : |x| = 1 or \y\ = 1} = $ _1 (0). (13) 

It follows from the computation of d§ / dt that the diagonal set A of T x T is repelling, and that the boundary 
cycle dT> is an attracting limit cycle. o 

The action of the vector field is to descend off of the "fins" of C onto the region T>, and then to circulate 
about while pushing out to the boundary cycle dT>, as in Figure |[ 

This example illustrates how one can use a relatively simple vector field on the configuration space to 
construct a pattern which is free from collisions. In fact, one could use this circulating flow to regularize 
potential collisions between AGV's in a general graph setting by localizing the dynamics near a pairwise 
collision to those on a trivalent subgraph. 
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4 Patterns and Vector Fields for Monotone Cycles 

Optimization of patterns in the workspace is deeply entwined with the geometry of the configuration space: 
in [ [Tot , it is shown that various Finsler structures on C can be chosen to optimize total distance traveled or 
net time elapsed. The net result of this inquiry is that minimizing Euclidean distance (in the product of 
the graph metric) on the cells of the configuration space yields locally optimal configuration sequences with 
respect to both distance traversed and elapsed time. In this section, we consider the problem of constructing 
vector fields which are tuned to trace out specific patterns of cooperative dynamics. We begin with a 
specification of a suitable language for describing patterns. 

4.1 A Grammar for Patterns 

The setting we envisage is as follows: the three ends of the graph T are stations at which an AGV can 
perform some function. The AGV pair is required to execute an ordered sequence of functions, requiring an 
interleaved sequence of visitations. In order to proceed with vector field controls for cooperative patterns, it 
is helpful to construct the appropriate symbolic language, as done in §|| for single AGV systems. Denote the 
pair of AGV states as x and y respectively. Also, denote the three docking stations as vertices v\ through 
V3 as in Figure |. The grammar Q we use is defined as follows: 

• (xi): These represent configurations for which the AGV x is docked at the vertex v i} i = 1..3. The 
AGV y is at an unspecified undocked position. 

• (yi): These represent configurations for which the AGV y is docked at the vertex Uj, i — 1..3. The 
AGV x is at an unspecified undocked position. 

• (xiyj): These represent configurations for which the AGV x is docked at vertex Uj while the AGV y is 
simultaneously docked at the vertex Vj, j =/= i. 

For example, the word (xl)(y2)(x3y2) executes a sequence in which the first AGV docks at Station v\, then 
undocks while the second AGV docks at Station v%. Finally, the AGV's simultaneously dock at Stations U3 
and z>2 respectively. 

As we have assumed from the beginning, the one-dimensional nature of the graph-constraints precludes 
the presence of multiple agents at a single docking station; hence, there are exactly twelve symbols in the 
grammar Q. From this assumption, it follows that particular attention is to be paid to those trajectories 
which do not make excursions onto the "fins" of the configuration space. It is obvious from the physical nature 
of the problem that planning paths which involve traveling on the fins is not a locally optimal trajectory 
with respect to minimizing distance or elapsed time. Suffice to say that we restrict attention for the moment 
to trajectories, and limit cycles for patterns in particular, which are constrained to the region T> C C. 

We identify each symbol with a region of the boundary of the unbranched portion of C; namely, dV is 
partitioned into twelve docking zones as in Figure |[ Note further that there is a cyclic ordering, <g, on Q 
induced by the orientation on the boundary of the disc along which the zones lie. By a cyclic ordering, we 
mean a way of determining whether a point q lies between any ordered pair of points (pi,P2)- 

We proceed with the analysis of limit cycles on C. Consider the class of pattern vector fields, Xp, on C: 
for every X G Xp, 

1. X defines a semiflow on C and a genuine flow off of the non-manifold set of C; 

2. There is a unique limit cycle 7 which is attracting and which traces out a nonempty word in the 
grammar Q; 

3. The diagonal set A is a repellor with respect to X; 

4. There are no fixed invariant sets of X which attract a subset of positive measure save 7. 

Then, the class of monotone vector fields, Xm, is that subset of Xp for which the limit cycle, 7, lies in T>. 
A word w composed of elements w = w\W2...w n in the grammar Q said to be monotone with respect to the 
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(x3y2) 




(x2y3) 

Figure 3: The six pairs of contiguous edges in dT> each corresponds to a configuration where one AGV is 
docked at an extreme vertex of the graph. Outermost vertices of dV are points where both AGV's are 
docked. Labelling the edges (xi) and (yj) and the vertices (xiyj) yieldsd the cyclically ordered grammar Q. 



cyclic ordering <g if Wi-\<Wi<Wi+i for every i (index operations all mod n). The following result justifies 
our use of the term monotone in describing those limit cycles which lie on the disc. 

Theorem 3 Within the class of vector fields Xm, the limit cycles trace out monotone words in the cyclically 
ordered grammar (Q, <a). 

Proof: Any limit cycle of the flow must be embedded (the curve does not intersect itself). After a 
small perturbation, one may assume that the boundary zone dT> is visited by 7 in a finite number of points, 
Q := jCidV. There is a cyclic order < t , defined via time with respect to the dynamics of the limit cycle: i.e., 
the order in which points are visited by 7. This is contrasted with the induced cyclic ordering <g on the set 
Q given by the orientation of the boundary curve dT>, up to a choice between clockwise and counterclockwise. 
The theorem follows from showing that <f = <a up to a cyclic permutation and a choice of orientation of 
dV. 

Induct on J the number of points in Q. For J = 1, the theorem is trivially true, so assume that the 
orderings are equivalent for all embedded curves on a disc with less than J boundary-intersections. Consider 
two points p,q 6 Q which are consecutive in the <t ordering. There is an embedded sub-arc a C 7 which 
connects ptoq within the interior of T>. By the Jordan Curve Theorem ||l7f , a separates T> into two topological 
discs; hence, 7 must lie entirely within the closure of one of these discs. Consider the open subdisc whose 
intersection with 7 is empty and collapse the closure of this disc to a point: this yields a modified curve 7 
which is still embedded in a disc having J — 1 intersections with the boundary, as illustrated in Figure ^. 
By induction, the ordering < t equals <q up to orientation on this subdisc. Reinserting the distinct points p 
and q by "blowing up" the crushed disc does not change the ordering properties, since these were chosen to 
be adjacent. o 

Hence, the ony admissible words in the grammar Q are those which are monotone. It is, however, possible 
to realize many if not all of the non-monotone cycles as limit cycles for a semiflow on the full configuration 
space C; one must design the semiflow so as to utilize the fins for "jumping" over regions of T> cut off by 
the limit cycle. Such vector fields quickly become very convoluted, even for relatively simple non-monotone 
limit cycles. 
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Figure 4: The embedded arc a divides T> (pictured as a smooth disc) into two discs, one of which is collapsed 
to a point x. 

4.2 Isotopy Classes of Limit Cycles 

Given a limit cycle 7 which traces out a pattern by visiting the boundary zone dT> in the ordered set Q C dT>, 
one wants to know which other limit cycles minimize a given performance functional while still visiting Q 
in the proper sequence. The mathematical framework for dealing with this problem is the notion of isotopy 
classes of curves. 

Two subsets A and Ax of a set B are said to be (ambiently) isotopic rel C (where C C B) if there exists 
a continuous 1-parameter family of homeomorphisms, f t : B~>B such that 

1. /o is the identity map on B\ 

2. f 1 (A Q )=A 1 ; and 

3. ft\ c is the identity map on C for all t. 

As £ increases, f t deforms B, pushing A a to A\ without cutting or tearing the spaces and without disturbing 
C. 

There are two ways in which optimization questions relate to isotopy classes of limit cycles: (1) Given 
an element of the grammar Q, in which isotopy class (rel the docking zones) of curves does an optimal limit 
cycle reside? (2) Within a given isotopy class of cycles rel Q, which particular cycle is optimal? 

For a monotone limit cycle on T>, the above question (1) focuses on the location of the cycle with respect 
to the central point (0,0), which is deleted from the disc T>. It is a standard fact from planar topology 
that every curve in the punctured disc has a well-defined winding number, which measures how many times 
the cycle goes about the origin, and, furthermore, that this number is either -1, 0, or 1 if the cycle is an 
embedded curve. This winding number determines the isotopy class of the curve in D. Hence, the problem 
presents itself: given an element of the grammar Q, which isotopy class rel the docking zones is optimal (with 
respect to any/all of the functionals defined)? Is the winding number zero or nonzero?^ 

To briefly address this question, we define the gap angles associated to a limit cycle. Given the docking 
zones Q — {qi, Q2, • • ■ , Qj} ordered with respect to time, we define the gap angles to be the successive 
differences in the angular coordinates of the qf. thus gaj := P{qj + \) — P(qj) 7 where P denotes projection of 
points in V onto their angular coordinates and subtraction is performed with respect to the orientation on 
dV. 

For simplicity, we consider the optimization— isotopy problem in the case of a discrete cost functional W<j, 
defined to be the intersection number of the path with the branch locus of C — i.e., the number of times an 
AGV occupies the central vertex (the shared resource in the problem). Similar arguments are possible for 
other natural performance metrics JTo| | . 

Proposition 5 There is a W ^-minimizing embedded monotone cycle on T> (rel a given docking zone Q) 
having winding number zero with respect to the origin if there is a gap angle greater than n. Conversely, if 
there are no gap angles greater than it, then there is a W ' ^-minimizing embedded cycle of index ±1. 

4 The difference between +1 and -1 is the orientation of time. 
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Proof: Let Q consist of the points {qj}{ on the boundary circle. The gap angles {ga,j}{ are the 
differences of the angles between the points qj and qj + i (indices mod J). Since J2j 9 a j = 2tt, there can be 
at most one gap angle greater than it. To simplify the problem, use a 1-parametcr family P t of maps from 
the identity P to the projection P = P 1 which deforms V to the boundary circle S := dV by projecting 
continuously along radial lines. The index of a curve on T> is invariant under this deformation, as is the 
function W^. 

Denote by 7^, the subarc of 7 between points qj and qj+i (all indices mod J). Denote by ctj the subarc 
of the boundary S between points qj and qj+i, where the arc is chosen to subtend the gap angle gcij. Since 
the boundary curve S = UjO? is a curve of index ±1, the arcs jj and ay are isotopic in T> rel their endpoints 
for all j if and only if 7 is a curve of index ±1. 

Assume first that there is a gap angle gaj > tt with 7 an index ±1 curve on S which intersects the branch 
angles = {nir/3 : n G Z} in a minimal number of points among all other closed curves on S which visit 
the points Q in the specified order. It follows that the arc P(jj) subtends an angle greater than w and thus 
increments by at least three. One may replace jj by a curve j'j which substitutes for the arc 7, one 
which wraps around 'the other way' monotonically. This changes the index of 7 from nonzero to zero, since 
the arc jj is no longer isotopic to ctj. Also, it is clear that this either decreases the number of intersections 
with 9 or leaves this number unchanged. 

We must show that the replacement arc 7^ can be chosen in such a way that is does not intersect the 
remainder of 7. However, since 7 is a curve of index ±1, we may isotope each arc 7, to the boundary curve 
cti without changing the value of W^. Thus, we may remove 7, and replace it with the curve which is, say, 
a geodesic (in the natural metric geometry) from qj to qj+i- As this curve does not approach the boundary 
S apart from its ends, the new curve 7' is an embedded curve of index zero without an increase in Wj. 

Assume now, on the contrary, that 7 is a minimizcr of index zero which has all gap angles strictly 
less than tt. Then each arc from 7, must intersect the branch set in at most three components, since, 
otherwise, the subtended arc would be in excess of 47r/3. In the case where there exists an arc with exactly 
three intersections with the branch set, this arc may be replaced by an arc which goes around the singularity 
in the other direction without changing the number of intersections with the branch set (since there are a 
total of six branch lines); however, the index of the curve is toggled between zero and nonzero. 

The final case is that in which each arc intersects the branch set in at most two places. However, since 7 
is a curve of index zero, some arc 7j must not be isotopic to ctj. Hence, the projection deformations P t must 
push 7j to a curve in the boundary which S whose subtended gap angle is 2ir — gaj > it. Thus, jj intersects 
the branch set in at least three places, yielding a contradiction. Replacing 7, by the appropriate arc which 
is isotopic to ctj yields a W^-minimal cycle of nonzero index. o 

4.3 Tuning Cycles 

In order to proceed with the construction of vector fields which realize monotone cycles, we work with vector 
fields on the smooth unit disc in R 2 and map these to the annular region T> of the configuration space via the 
push-forward induced by the natural homeomorphism. It will be convenient to keep track of which "wedge" 
of the annular region a point (r, 9) is. To do so, we introduce a parity function 

P(9) := (-1)U30/tJ+|6<?AJ} 5 ( 14 ) 

where [t\ is the integer-valued floor function. Recall the notation for the intrinsic coordinates for a point x 
on the graph T: x — \x\e L ^ x ^ where \x\ G [0, 1] is the distance from x to the central vertex, and e t ( x ) is the 
unit tangent vector pointing along the direction of the t(x)-edge. Here the index, l(x) is an integer (defined 
modulo 3) and will be undefined in the case when \x\ = 0, i.e., x is at the central vertex. 

Lemma 6 The following is a piecewise-linear homeomorphism from the punctured unit disc in R 2 to the 
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subset T>: F(r, 9) = (x, y) where, 



\v\ = 

l(x) 

<y) 



cot %d 



Z7T 



V{9) 

vie) 



(15) 



T/ie inverse of this homeomorphism is given by F (x,y) = (r,9), where, 




r = 



3 tan 

V{9) = +1 



0-i) 



■'(I/) + 1 



= t(at) + 1 
or | re | = 
l(x) = 



\x\ 



or \y\ = 



(16) 



Note that all 9 values are defined modulo 2ir and all index values are integers defined modulo 3. 

Proof: Begin by working on the region "Dtp C T> where i(x) = 1 and i(y) = 2. As noted earlier, this 
subspace is isometric to the positive unit square in R 2 with the origin removed. We need to map this to 
the subset {(r, 9) : r £ (0, 1], 9 € [0, 7r/3]}. The simplest such homeomorphism is to first shrink along radial 
lines, leaving the angle invariant; hence 



F 

\u\ 



\x\ < \y\ 
\y\ < \x\ 



(17) 



Next, we squeeze the quarter-circle into a sixth of a circle by multiplying the angle by 2/3, leaving the radial 
coordinate invariant: 

^Itan^M. (18) 
3 \x\ 

This gives the basic form of F~ x as per Equation (|l6|) . To extend this to the remainder of V, it is necessary to 
carefully keep track of l(x) and t(y) and subtract off the appropriate angle from the computation of 9. Also, 
the condition of |x| < \y\, etc., in Equation ( |l7j ) is incorrect on other domains of T>, since the inequalities 
flip as one traverses from square to square: the parity function V{9) keeps track of which "wedge" one is 
working on. 

To determine F from F~ x is a tedious but unenlightening calculation, made more unpleasant by the 
various indices to be kept track of. Briefly, given r and 9 on the first sixth of the unit disc, one knows from 
Equation ( [l7|) that either \x\ — r or \y\ — r, depending on whether 9 is above or below 7r/4. To solve for 

respectively. To 



tanH6> 



or \x\ 



cot : 



the other magnitude, one inverts Equation (18) to obtain \y\ = r 
generalize this to the other Vij domains of T>, it is necessary to take absolute values and to use the parity 
function V(9) as before. Finally, the computation of the index is obtainable from the combinatorics of the 
coordinate system as illustrated in Figure @. o 



Hence, by taking the push-forward of a vector field X = (f, 
smooth vector field, 




tan(|0) 
cot ( | (9) 



+ |r0 scc 2 (|e») 
+ §r0 csc 2 (§6>) 



with respect to F, one obtains the piecewise 
V{9)=+1 
) V(9) = -l 



(19) 
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which simplifies to: 



/ M 
\y\ = ' 

V 

( lil 
V \y\ 



.11 


- + ¥— 


\x\ 


X 


14 


(\y\V 






[mj 


X 


- + ¥— 


\y\ 


y 




/i i\ a 

/ \x\ \ 



V(6) 



V(d) 



(20) 



Given a simple closed curve 7 in R 2 which has nonzero winding number with respect to the origin, 7 
may be parametrized as {(r, 9) : r — /(#)} for some periodic positive function /. To construct a vector 
field on R 2 whose limit sets consist of the origin as a source and 7 as an attracting limit cycle, it suffices 
to take the push-forward of the vector field r — r(l — r) 9 = lo under the planar homeomorphism 
4> : (r,9) 1— > (f(9)r,9), which rescales linearly in the angular component. The calculations follow: 



D<t> 



r 1 - 



. r 
f 

r-fio 



f rf> 
1 



r(l - r) 
10 



r 

7 



(21) 



Hence, given f(9), we may tune a vector field to trace out the desired limit cycle and then use Equations ( |l5| ) 
and (|l^) to map it into intrinsic coordinates. 



4.4 Optimal Chords within a Hybrid Controller 

To design optimal cycles with winding number zero, then, we turn to constructing customized portions of 
limit cycles, or chords which can be pieced together via a state-actuated hybrid controller, much as in §^|. In 
other words, instead of building a simple fixed vector field with a limit cycle, we will use a set of vector fields 
which vary discretely in time and which may be pieced together so as to tune a limit cycle to the desired 
specifications. There is nothing in this construction which relies on the index-zero property and thus these 
chords can be used to generate all monotone limit cycles on C. 

Let G denote a word representing a desired monotone limit cycle on the configurations space C. Choose 
points {qi} on the boundary of T> which correspond to the docking zones for the cycle given by G. Choose 
arcs Ui on T> which connect qi to qi+\ (using cyclic index notation). The arcs Ui are assumed given in the 
intrinsic coordinates on T>, as would be the case if one were determining a length-minimizing curve. 

In the case where the limit cycle a :— UiCti is an embedded curve of nonzero index, the procedure of the 
previous subsection determines a vector field X a on C which realizes a as an attracting limit cycle with the 
appropriate dynamics on the complementary region. Recall: one translates a to a curve on the disc model 
via the homeomorphism of Equation ([l6|). Then, representing the limit cycle a as a function f a (9), one takes 
the vector field of Equation ( |2l| ) and, if desired, takes the image of this vector field under Equation (p(i|). 

If, however, this is not the case, consider the arc ctj for a fixed j and construct an index ±1 cycle 
Z?- 7 = Uif3f which has docking zones {qi} such that /3j = ctj. Then the vector field as constructed above 
has (3 as an attracting limit cycle. Denote by <i> J the Lyapunov function which measures proximity to (3: 
:— — (with distance measured in say the product metric on C). Then, consider the modified 
Lyapunov function $ J '(p) := &(p) + \\p — qj+i\\ , which measures the distance to the endpoint of the arc /3j 
in addition to the proximity to /3 3 . 

Repeat this procedure for each j, yielding the vector fields {X^} which attract respectively to limit cycles 
{P . It follows that X^ prepares X^ +1 since the goal point of X^ , qj + i lies on the attracting set of X^ +1 . 
The Lyapunov functions {>& J } serve as a set of funnels which channel the orbit into the sequence of arcs ctj, 
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forming a. One scales the ty 3 so that a ^ < e event triggers the switching in the hybrid controller from X 3 



By construction, the hybrid controller ( p2| ) realizes a limit cycle within e of a as the attracting set. 

5 Future Directions 

A point of primary concern is the adaptability of the global topological approach to systems which increase 
in complexity, either through more intricate graphs or through increased numbers of AGV's. The latter 
is of greater difficulty than the former, since the dimension of the resulting configuration space is equal to 
the number of AGV's. Hence, no matter how simple the underlying graph is, a system with ten indepen- 
dent AGV's will require a dynamical controller on a (topologically complicated) ten-dimensional space: a 
formidable problem both from the topological, dynamical, and computational viewpoints. 

However, there are some approaches which may facilitate working with such spaces. Consider the model 
space C with which this paper is concerned: although a two-dimensional space, C can be realized as the 
product of a graph (a circle with six radial edges attached) with the interval (0, 1]. In fact, if we consider the 
circulating flow of Equations (@)-(^), one can view this as a product field of a semiflow on the graph (which 
"circulates") with a vector field on the factor (0, 1] (which "pushes out" to the boundary). 

A similar approach is feasible for arbitrary graphs. The following result has recently been proven S: 

Theorem 4 Given a graph T, the configuration space of N distinct points on T can be deformation retracted 
to a subcomplex whose dimension is bounded above by the number of vertices of T of valency greater than 



This theorem implies the existence of low-dimensional spines which carry all of the topology of the 
configuration space. For example, the above theorem implies that the configuration space of N points on 
the Y-graph can be continuously deformed to a one-dimensional graph, regardless of the size of N. Since 
the full space can be deformation retracted onto the spine, a vector field defined on the spine can be pulled 
back continuously to the full configuration space, thus opening up the possibility of reducing the control 
problem to that on a simpler space. Addtional results about the topology of configuration spaces on graphs 
may yield computationally tractible means of dealing with complex path planning: for example, having a 
presentation for the fundamental group of a configuration space of a graph in terms of a suitably simple 
set of cycles would be extremely well-suited to a hybrid control algorithm based on "localized" vector fields 
supported on small portions of the full configuration space. 

The optimization problem is another avenue for inquiry. The fact that a dynamical approach allows 
for increased density of AGV's on a graph (as compared with blocking-zone strategies) would indicate an 
increased efficiency with respect to, say, elapsed time-of- flight. However, a more careful investigation of the 
tuning of optimal cycles is warranted. 

We believe that the benefits associated with using the full configuration space to tune optimal dynamical 
cycles justifies a careful exploration of these challenging spaces. 

A The Topology and Dynamics of Graphs 

In this appendix, we provide a careful basis for the use of vector fields on configuration spaces of graphs. 
In the setting of manifolds, all of the constructions used in this paper are entirely natural and well-defined. 
However, on spaces like C, the most fundamental of notions (like the Existence and Uniqueness Theorems 
for ODEs) are not in general valid. 

We begin by defining vector fields on graphs. For present purposes, we find it convenient to work with 
an intrinsic formulation (i.e., directly in the graph rather than via an embedding in a Euclidean space) of 

5 We have since learned of two others who have independently proved this result: fujl, having learned of these spaces fron 
[nil; and 111, who discovered these spaces while working on the topology of Brownian motion on graphs. 



to XJ+ 1 : 




(22) 
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these objects. To this end, denote by v a vertex with K incident edges {ei}f , and by {Aijf- a collection 
of nonsingular vector fields locally defined on a neighborhood of the endpoint of each a (homeomorphic to 
[0,1)). 

Lemma 7 A set of nonsingular vector fields {Xi} on the local edge set of a graph T generates a well-defined 
semiflow on T if 

1. Each edge field Xi generates a well-defined local semiflow on (0,1); and 

2. The magnitude of the endpoint vectors ||Xi(0)|| (taken with respect to the attaching homeomorphisms) 
are all identical; and 

3. Among the signs of the endpoint vectors Xi(0) (either positive if pointing into [0, e) or negative if 
pointing out) there is a single positive sign. 

Proof: Since the vector field is well-defined away from the vertex, it is only necessary to have the 
magnitudes ||Xi(0)|| agree in order to have a well-defined function ||X|| on T. In order to make this a 
well-defined field of directions, we must also consider in which direction the vector is pointing. Again, this 
is determined off of the vertex by (1). Condition (3) means that at the vertex, there is a unique direction 
along which the vector field is pointing out: all other edges point in. Hence, the direction field, as well as 
the magnitude field, is well-defined. 

The semiflow property follows naturally from this. Assume that the N th edge of T has the positive sign. 
Then, given an initial point x S T, if x S ejv, then the orbit of x under the local field Xjq remains in e^v and 
is well-defined. If x G ej for some j ^ N, then the union of the edges U ejv is a manifold homeomorphic 
to R on which the vector fields Xj and Xn combine to yield a well-defined vector field, since the directions 
are "opposite." As we are now on a manifold, the standard Existence Theorem implies that x has a forward 
orbit (which passes through the vertex and continues into e^r). Thus every point on T has a well-defined 
forward orbit. o 

In the case where the vector fields have singularities, it is a simpler matter. If the singularities are not at 
the vertex, then there is no difference. If there is a singularity at the vertex, then condition (3) in Lemma |?] 
is void — all such vector fields are well-defined. 

In order to extend these results to the configuration space of this paper, consider the space C — T x T — A 
and let (x, y) G C denote a point on the branch set of C. Because of the structure of T and the fact that the 
diagonal points are deleted, it follows that at most one AGV may occupy a non-manifold point of T. Hence, 
a neighborhood of (x, y) in C has a natural product structure N = T x R. Let P : N — > T denote projection 
onto the first factor. 

Lemma 8 A nonsingular vector field X on the individual cells of C generates a well-defined semiflow if (1) 
the projection of the local vector fields onto the graph factor, P*(X\^ xT ), satisfies Lemma^for each point 
x in the branch set of C; and (2) the projections of the vector fields on the branch set to the R-f actor are 
equal up to the attaching maps. 

Proof: Off of the branch set, the space is a manifold and hence the vector field gives a well-defined flow. 
If p is a point on the branch line, condition (2) implieds that the vector field is well-defined with respect 
to the attaching maps and the net effect in the i?-factor is a drift in this direction. In the graph factor, 
condition (1) and the proof of Lemma implies that there is a unique forward orbit through p. o 

Heuristically, this condition means that, as in the case of a graph, the vector field must point "in" on all 
but one sheet of the configuration space in order to have well-defined orbits. We may thus lift the criteria 
of Lemma [j] to the product configuration space. All of the vector fields in this paper are so constructed. 
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